Spin contribution to the ponderomotive force in a plasma 



O 



G. Brodin, A. P. Misra, and M. Markluncfl 

Department of Physics, Umed University, SE-901 87 Umea, Sweden 
(Dated: 24 March, 2010; Revised 27 July, 2010) 

The concept of a ponderomotive force due to the intrinsic spin of electrons is developed. An 
expression containing both the classical as well as the spin-induced ponderomotive force is derived. 
The results are used to demonstrate that an electromagnetic pulse can induce a spin-polarized 
plasma. Furthermore, it is shown that for certain parameters, the nonlinear back-reaction on the 
electromagnetic pulse from the spin magnetization current can be larger than that from the classical 
free current. Suitable parameter values for a direct test of this effect are presented. 
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The use of the spin properties of material constituents 
for e.g., carrying information is currently an important 
paradigm |1J]. However, the spin properties of the ma- 
terial constituents also make its presence felt through 
collective effects. In particular, recent findings point to 
the possibility of observing quantum plasma effects [2] 
through the electron spin [3] in regimes otherwise thought 
to be classical [4] . Such results are due to the complex in- 
terplay between collective plasma effects and the system 
nonlinearities. In classical plasmas, nonlinear effects play 
an important, sometimes a crucial, role. For example, 
the density fluctuations induced by the ponderomotive 
force of an electromagnetic (EM) wave lead to an elec- 
trostatic wake field [5] , as used in advanced particle accel- 
erator schemes [6]. In other regimes, the back-reaction 
on the EM-wave due to the density fluctuations leads 
to phenomena such as soliton formation, self-focusing or 
wave collapse 0. Such radiation pressure-like effects 
are widely used in high-intensity laser experiments Q, 
and generalizations to include certain types of quantum 
plasma effects have recently been made However, to 
our knowledge the possibility of spin induced contribu- 
tion to the ponderomotive forces has not been explored. 

In the present work we will solve the full set of equa- 
tions for the spin dynamics of charged particles in the 
presence of a weakly nonlinear EM wave pulse, propa- 
gating parallel to an external magnetic field, in order to 
find the contribution to the ponderomotive force. In the 
classical limit, we recover the well-known expression first 
derived by Karpman and Washimi The spin con- 

tribution to the ponderomotive force will in general act 
in opposite directions for spin-up and spin-down popu- 
lations. As a consequence, an EM-pulse (due to, e.g., a 
laser or a microwave source) may induce a spin-polarized 
plasma. In particular, it is demonstrated that this mech- 
anism can induce large spin-polarization for a laser source 
in the UV-regime. For this case it should be noted that 
the effect of the external magnetic field is negligible as 
the laser frequency is much higher than the cyclotron fre- 
quency, but our general expression applies also for low- 
frequency (If) waves in magnetized plasmas. 

When combined with the high-frequency (hf) oscilla- 



tions, the classical If density response generates a hf cur- 
rent that results in a cubical nonlinearity. This classical 
nonlinearity can be compared to a cubically nonlinear 
magnetization current, caused by the electrons spins in 
a plasma that are spin-polarized due to ponderomotive 
effects. It turns out that for a plasma frequency cor- 
responding to a metal density and a hf source (e.g., an 
x-ray free electron laser (XFEL)), the spin contribution 
can be larger than the classical contribution. 

We will assume the existence of two electron popula- 
tions, namely spin-up and spin-down relative to a back- 
ground magnetic field Bp = Sqz, to be denoted by u and 
d respectively, and formally treated as different species. 
Such division are relatively common in e.g. semiconduc- 
tor physics Although the spin states of the parti- 
cles will be perturbed by the presence of electromagnetic 
waves, the separation of species is still well-defined pro- 
vided the physics associated with spin-flips can be ne- 
glected (see Refs. Q and 11 1 for a discussion and a cal- 



culation of spin- flip frequencies, respectively). The basic 
equations take the form 0, |3| 

dtUa + V • (UqVq) = 0, (1) 

m {dt -I- Vq • V) Vq = g (E + Vq X B) - WPa/ria 

+ {^^i/h)SlVBa. (2) 
{dt + V, • V) = -(2M/ft)B X S„, (3) 

and V • E = (q/eo) (?^u + — rii), where Sq is the spin 
of species a (with a = u,d), rii is the ion density, and 
q = — e is the electron charge, /i = —9/13/2, where 
fiB = eh/2m is the Bohr magneton, g « 2.0023192 is the 
electron (/-factor, and Pa is the pressure [l2|. The princi- 
pal condition for the validity of Eqs. (UJ-Q is that char- 
acteristic spatial scale lengths is longer than the thermal 
de Broglie length. Naturally Eqs. ([I])-© also has the in- 
herent limitations of a fluid model, i.e. resonance effects 
for particle velocities matching the phase velocity and/or 
group velocity are left out. Furthermore, effects due to 
the off-diagonal componenets of the pressure tensor are 
omitted. In what follows we consider an electron-proton 
plasma with classical ion dynamics. 

Our derivation of the ponderomotive force below will 
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be based on a direct perturbative analysis of Eqs. ([T])- 
([3]). We note that _an alternative approach based on 



Lagrangian method 1^ [16| can be useful. Assuming a 
slowly varying plane EM wave E =Eexp[i{kz—Lut)]+c.c., 
where c.c. denotes complex conjugate, Eq. ^ becomes 



{q/ni)'Ei — rivc X z, 



(4) 



where = eBo/m is the electron-cyclotron frequency. 
Here we have assumed that the background state con- 
tains no net drifts. Note that kinetic effects, which is 
of potential importance (see e.g. [Tt'I) but is outside the 
present model, is assumed to be negligible throughout the 
manuscript. Within a fluid model, the average longitudi- 
nal velocity, that mOT contribute through the convective 
derivative (cf. Ref. [iSj ) becomes second order in the am- 
plitude, and hence can be neglected here. For notational 
convenience, we have dropped the tilde denoting the en- 
velope function, and it is understood that all derivatives 
act on the slowly varying amplitudes. Next, we define 



the variables Va± 



± ivay, E± = ± iEy. Sub- 



stituting the lowest order result Va± = iqE±/m (w ± fl) 
into the correction term in Eq. we obtain 



V± = Va± 



iE± 



1 dE± 



uj±n dt 



(5) 



Using Faraday's law, V x E = — 9tB, we similarly obtain 
the expression for the perturbed magnetic field as 



dz 



dt 



(6) 



The classical ponderomotive force component is 



F,, 



= ^(vxB), = 
m 




for RCP, 



for LCP. 



Substitution of Eqs. (O and (O into Eq. ([7]) gives 



d 



± 



d 



2vn?u> {uj ±0) \_dz w (cj ± 17) 9i 



\E\ 



(7) 



(8) 



in agreement with the classical result |14l |. 

Next, we derive the effects due to the finite magnetic 
moment of the electrons. Through the force Faz = 
(2fi/mh){S^VBa)z in the averaged momentum equation, 
a ponderomotive effect due to spin will be generated, 
where the EM wave will lead to the separation of the 
spin-up and down electrons, as will be shown below. 
Starting from the linearized spin-evolution equation 

{dt - zw) S„ = -{2^i/n) {Boz xS^ + SoaBxz), (9) 

where Squ — ti/2 — —Sq^, the contribution from the 
magnetic dipole force can be obtained. Neglecting the 
slow time derivative, Eq. © gives 

Sa± = Sax ± iSay = T'^^lSQaB± / [H (w ± UJg)], (10) 



where u)g = gfiBBo/h = {g/2)Q is the spin-precession 
frequency. Then, including the first order correction, the 
expression for the perturbed spin becomes 



Sa± — 



2/iS'oo 



h{uj ± UJg) 



tb±± 



dB± 



(w ± cjg ) dt 



The spin-ponderomotive force can be written as 

Fa, = {2^i/mh){Sa±dzBl + Sl±d,B±). 
Substitution of Eq. ([11]) into Eq. ([12]) gives 



4^^ Spa 

mh? (cj ± Wg) 



d 



d 



dz (uj ±ujg) dt 



\B\ 



(11) 



(12) 



(13) 



The above expression applies to arbitrary EM wave prop- 
agation parallel to Bq. The overall structure of the spin 
ponderomotive part of the force is similar to the 
classical part However, there are differences. Firstly, 
the frequency resonances occur at the spin precession fre- 
quency UJg = {g/2)n 1.0011617. Secondly, the depen- 
dence on the unperturbed spin state means that spin-up 
and spin-down populations are pushed in opposite direc- 
tions by the spin force. Thirdly, for frequencies well be- 
low the cyclotron frequency typically the part of the spin 
contribution proportional to the time-derivative is negli- 
gible, whereas it is crucial for the classical contribution. 
For frequencies well above the cyclotron frequency, the 
force ratio scaling is \Faz\ / \Fcz\ = hk{l + Vg / Vp) /mvp ~ 
hk/mc for Vg,Vp c, where Vg(^p) is the group (phase) 
speed of the wave. As we will see below, even a rather 
weak spin-ponderomotive force, corresponding to moder- 
ately high frequencies, can lead to large modifications of 
the nonlinear dynamics in an unmagnetized plasma. 

We now use the expressions for the ponderomotive 
forces as source terms for longitudinal If perturbations. 
We define iVi.2 = riu ± rid and Vi^2 = {vu ± Vd) /2. In 
what follows, we will also neglect any difference in the 
unperturbed spin populations, i.e., we will use no„ = 
nod = no/2, which is a good approximation when the 
Zeeman energy is smaller than the thermal energy. From 
the If parts of the continuity equations for spin-up (u) 
and spin-down (d) populations we then have 



dtN, 



-npdzVi,2, 



(14) 



From the momentum balance equations, using Eq. ([6]), 
we obtain for the If response the equations 



dVi ^ 

dt m 2m?uj {uj ± il) 



d 



± 



kn 



d 



dz UJ {uj ± Q) dt 



\E[ 



and 

dt 



Vfc^S-o 



mh?uj'^ (uj ± UJg) 



d 



d 



dz {uj ± UJg) dt 



\E' 



(15) 
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where Ei is the If part of the electric field, Sq = h/2 
and ujp = (nog^/meo)^^^- Here we have neglected ther- 
mal effects, which is justified if v^y^ <^ v^, where Wth is 
the thermal velocity. Similarly particle dispersive effects 
and Fermi pressure effects, which may influence the lon- 
gitudinal dynamics has been neglected, which can 
be justified in the example considered below [l^. With 
immobile positive charge carriers, Poisson's equation is 
d,Ei = (g/eo)(iVi - no). Together with Eqs. (fT4 | - (fT5|) . 
we then obtain the wave equations 



2muj {uj ± Vl) 



It 



kVgil 



d^\El_ 

(16) 



= ± 



kvn 



iuJ±UJg) 



d^\E[ 



(17) 

where we have transformed to a comoving frame, with 
^ = z — Vgt. The spin polarization [see Eq. pT)) ]. can 
be integrated directly to give N2 oc whereas A^i 

is non- locally related to \E\'^ through (fT6| . due to the 
possible excitation of a plasma oscillation wakefield with 
a characteristic wavelength Xp = Vg/ujp. To demon- 
strate that the spin effects can be significant also when 
Bq 0, we compare the amplitude of the total den- 
sity perturbation A^i with the degree of spin-polarization 



N2 in an unmagnetized plasma ujg, il — 0. Further- 
more, to be specific, we consider hf EM waves with 
Vg,Vp ^ c. Finally, we use ujp < kc and use the es- 



timate d^\E\^/de --\E\^/Ll, where Lp > k'^ is the 
length of the hf pulse. The degree of spin-polarization 
is then N2/N1 ^ hijjp{kLp)'^ /mc^ . We note that the 
omission of ion density dynamics constrains this expres- 
sion to pulse lengths fulfilling Lp < c/ujpi, where ujpi 
is the ion plasma frequency. Thus, here we will con- 
sider the case of an EM-pulse interacting with a plasma 
without positive mobile charge carriers, i.e., a metal with 
uip/2TT ~ 10^^ s^^. A numerical example with a UV-laser 
of wavelength, A = 80 nm and pulse length, Lp = 15 /im 
leading to moderate spin-polarization ( A''2 /Ni « 3 at the 
centre) is displayed in Fig. 1. A longer pulse length or 
a shorter wavelength will give a higher degree of spin- 
polarization, i.e. a strongly spin-polarized plasma with 
N2 ^ A^i can be reached. We note that the polariza- 
tion of the EM wave is crucial. In the limit considered 
here (w 2> f2), the spin contribution to the ponderomo- 
tive force has opposite direction for RCP and LCP waves. 
Thus, an experiment along these lines must use circular 
polarized rather than linearly polarized light, as the effect 
is a factor il/uj smaller in the latter case. A related use of 
the ponderomotive force for isotope separation has been 
suggested in Ref . , where the different charge-to-mass 
ratios of different isotopes was used. 

Next, we compare the back-reaction on the EM-pulse, 
induced by the classical density perturbation A^i and its 



spin-polarized counterpart N2. From the classical current 
J = q{nuVu + UdVd) = qNiV, we have 

J± = qNiv± = iq^NiE±/[m {uj ± fl)], (18) 

with v„ — Vd for RCP and LCP waves. The 
magnetization current Jm — V x (M„ + 
Md) = (^/n) [V X (n„S„ +ndSd)], i.e., Jm± 
±{kgfj,B/2h){nuSu±+ndSd±), together with Eq. pU)) 
and the lowest order expression of B± [see Eq. ([6])], gives 

Jm± = ±&i (kfi/ujhf [ljSo/{uj ± Ug)] {N2E±) . (19) 

Now, Eq. (|17p can be integrated for N2 ■ A Gaussian 
pulse \E\ = £'0 exp(-^^/Lp) gives (for B^^LOg -> 0) 



m^Lj^Vg 



Furthermore, for a pulse length much larger than the 
plasma oscillation wavelength Ap, we can use the follow- 
ing estimate from Eq. as {Bq, Q ^ 0) 

N, ^ [eoOjl\Eo\^ eM-'^eiLl)]l{mu\lLlkl). (21) 

The density ratio for RCP and LCP waves is then given 
by 



N2 



\mc^ J \Vp 



k'^Vg 



1 + ^ 



(22) 

where kp = 1/Xp. The ratio of the two currents for 
RCP and LCP waves is then given by P = \Jm±/J±\ ~ 
{huj/mvl)\N2/Ni\, i.e. 

r (hojp/mc^y {kLp f [c^/vpVgY (1 + Vg/vp) . (23) 



Thus, for Vg,Vp ^ c, we have 



{fvjjp/mc^Y (kLp 



(24) 



Solving Eqs. ^ and (H?]), we in Fig. 2 compare the 
two current profiles for an XFEL with A = 1 nm, a pulse 
length Lp = 30 /zm, and a metallic plasma density, giv- 
ing ujp/2TT = 10^^ s~^. Our estimate ([M)) is then verified, 
and it is found that the central value of the current ratio 
is r ~ 3. These parameters are relevant for the XFEL 
at DESY 2l|. In fact, the shortest wavelength gener- 



ated by this facility is A = O.lnm, making the quantum 
mechanical back-reaction much larger than the classical 
response [P ~ 200 according to Eq. ([M]) ]. 

In the present Letter, we have generalized the classi- 
cal expression for the ponderomotive force in a magne- 
tized plasma to include the effect of the electron spin. 
Our main result, Eq. ()13p . applies for arbitrary electro- 
magnetic waves propagating along an external magnetic 
field. One of the main features of the spin-ponderomotive 
force is that it can induce a strong spin-polarization in a 
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FIG. 1. The normalized classical density perturbation A'^i 
and the spin-induced density difference N2 together with a 
Gaussian EM-pulse, \E\ = |i5o| exp(— ^^/Lp), as calculated 
from Eqs. (|16p and (|17p . The parameters correspond to an 
unmagnetized plasma with ijjp/2-K = 10^" s^^, A = 80 nm, and 
Lp = 15 /im (a.u.). 




FIG. 2. The profiles of the normalized spin-induced cur- 
rent density Jm and the classical current density J for a 
Gaussian EM-pulse, as obtained from Eqs. (|18l) and (|19|) . 
The parameters correspond to an unmagnetized plasma with 



U]p/2-K ■ 



A = 1 nm, and Lr, 



30 /im (a.u. 



plasma, even if the initial up- and down- states of elec- 
trons are equally populated. An example with an EM- 
pulse in the UV-regime is given in Fig. 1. Furthermore, 
even in an unmagnetized plasma, the nonlinear back- 
reaction from the spin-induced current can be larger than 
the classical back-reaction, provided the EM-pulse has a 
sufficiently short-wavelength. It should be stressed that 
we have here only compared the spin-polarized current 
response with the classical current density, and that rela- 
tivistic nonlinearities (see e.g. 0) that has been omitted 
here may play a role for the full dynamic evolution. An 
example with an XFEL is given in Fig. 2. Finally, we 



want to point out that the possibilities of nonlinear spin 
effects is still a relatively unexplored area, and general- 
izations, such as arbitrary directions of propagation, are 
likely to lead to new and interesting discoveries. 
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